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PROCESS FOR SIMULATION OF THE RESPONSE OF A DETECTOR OF 
RADIATION EMITTED BY RADIOACTIVE OBJECTS AND PROCESS 
FOR INSPECTION OF NUCLEAR FUEL ELEMENTS USING THIS 

SIMULATION 

DESCRIPTION 

Technical field 

This invention relates to a process for simulation 
of the response of a detector of radiation emitted by- 
radioactive objects and a process for inspection of 
5 nuclear fuel elements using this simulation. 

It is particularly applicable to the inspection of 
nuclear fuel rods that contain stacks of pellets of 
this fuel, these pellets emitting 7 radiation. 

10 State of prior art 

It is known how to inspect a set of rods of this 
type to verify the homogeneity of pellets in a given 
category . 

In order to do this, a preliminary calibration is 
15 made of the y radiation measurement system used for the 
inspection. For example, the detector for this system 
is described, in documents FR-A-2437002 , EP-A-0009450 
and JP-A-1527161 that includes an annular scintillator 
that is preferably made of sodium iodide activated with 
20 thallium Nal(Tl). 

This calibration also involves making specific 
rods and passing them in front of the detector, to 
obtain what are referred to in statistics as regression 
straight lines, and which in this case give the 
25 response of the detector to powder mixes used in the 
pellets depending on the content of their components 
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(for example uranium and plutonium) , for homogenous 
portions of . rods or for isolated pellets in one 
category among a group of pellets of another category. 

This requires a large number of measurements and 
5 an equally large number of calibration rods 
representative of these typical situations, but 
unusable in a reactor. 

Presentation of the invention 

10 The purpose of this invention is to overcome the 

disadvantages mentioned above by proposing a process 
that can economically control nuclear fuel rods, or 
more generally nuclear fuel elements. In this process, 
the preliminary calibration of the measurement system 

15 is eliminated and is replaced by a simulation of the 
response of the measurement system detector, in other 
words the count made by this measurement system. 

Specifically, the initial purpose of this 
invention is a process to simulate the response of a 

2 0 radiation detector detecting radiation emitted by 

radioactive objects, these objects containing 
radioelements or mixes of radioelements , this process 
being characterized in that: 

radioactive emission spectra representative of 
25 radioelements or mixes of radioelements are memorized, 

the detection characteristics of the detector are 
determined, 

the operating characteristics of received 
radiation are determined, 

3 0 - radioelements or mixes of radioelements among 

those for which spectra have been memorized are 
selected, and 
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the detection characteristics and operating 
characteristics are processed so as to individually 
reproduce radiation emitted for the chosen 
radioelements or mixes of radioelements, to obtain the 
simulated detector response. 

Preferably, the detection characteristics comprise 
data representative of the thickness through which the 
radiation passes before it is detected. 

Also preferably, operating characteristics also 
include the detector aperture 'angle, detected energy 
bands and electronic amplification characteristics of 
the detector. 

According to a preferred embodiment of the process 
according to the invention, regression straight lines 
are also built up starting from the simulated response. 

For example, the detector may be a y radiation 
detector and the invention is applicable particularly 
in the case in which the said objects are nuclear fuel 
elements . . 

The invention also relates to a process for the 
inspection of a set of nuclear fuel elements using the 
simulation process (applied to objects composed of 
these elements) , inspection process in which: 

the real composition of any of the elements of the 
assembly is analyzed, 

the detector is calibrated with this element for 
which the real composition has been analyzed, 

the simulated response is corrected using the 
response of the detector obtained during calibration, 
and 

all elements are inspected. 



B13067 PV 



For example, the elements may be nuclear fuel 
rods, these rods including stacks of pellets from this 
nuclear fuel . 

For example, in this case the detector may 
comprise an annular scintillator, and for example a 
sodium iodide scintillator may be used. 

Brief description of the drawings 

This invention will be better understood after 
reading the description of the example embodiments 
given below, which is for information only and is in no 
way restrictive, with reference to the attached 
drawings in which: 

FIGURE 1 is a schematic view of a detector for 
which the response is to be simulated, on a plane 
perpendicular to the axis of this detector, and 

FIGURE 2 is a schematic view of the detector shown 
in its shielded containment along a plane parallel to 
its axis. 

Detailed presentation of particular embodiments 

It is assumed that it is required to inspect a 
batch (in other words a set) of nuclear fuel rods. 
Each rod is a stack of pellets, for example containing 
uranium oxide and/or plutonium oxide. The pellets in 
these fuel rods are inspected individually. For 
example, this can be done using the detector described 
in the documents mentioned above. The structure of 
this detector is described with reference to FIGURES 1 
and 2 throughout the rest of . this description. 

It is a y radiation detector that comprises an 
annular-shaped scintillator 1 associated with three 
photomultipliers 2, 3 and 4. FIGURES 1 and 2 show a 
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rod .16 to be inspected that is composed o£ pellets 5. 
The detector also comprises a diaphragm or collimator 
6. This diaphragm limits the y radiation flux emitted 
by this pellet towards the scintillator, to 
5 approximately the length of each pellet. The three 
photomultipliers are uniformly distributed around the 
periphery of the scintillator. The outputs from these 
photomultipliers are connected to electronic 
measurement means forming a counting system 7, which 

10 will be described later. 

The scintillator is divided into identical sectors 
10, 11 and 12 that are optically isolated from each 
other and associated with photomultipliers 2, 3 and 4 
respectively. This scintillator is preferably of the 

15 sodium iodide type, activated with thallium. Layers 
13, 14 and 15 of an optical insulator such as aluminium 
can be seen, that optically isolate the different 
sectors in the scintillator. 

Figure 2 shows a shielded containment E protecting 

20 the detector against external y radiation that could 
disturb the measurements. Ifalso shows that pellets 5 
in rod 16 are contained in cladding 17 . This rod is 
moved along a direction 18 by means not shown. It also 
shows two annular parts 19 and 2 0 that form the 

2 5 diaphragm 6 and that are opaque to y radiation. The 
spacing between these parts is e, that can be adjusted 
by means not shown. < The rod to be inspected is moved 
along the axis 21 of the detector. 

The count system 7 comprises amplifiers- 

30 stabilizers 22 associated with photomultipliers, an 
adder 23 with inputs connected to these amplifiers, and 
single channel analyzers 24 (four in the example shown 
in figure 2) with inputs connected to the output from 
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adder 23, respectively. A computer 25 is provided for 
processing of signals output by these single channel 
analyzers 24. This computer is used with a memory 26 
and display means 27 and is also designed to run the 
5 software used to simulate the detector response. 

When a measurement is made on a pellet 5, a total 
spectrum is obtained for all isotopes in this pellet by 
summation (number of pulses per energy channel as a 
function of the energy) and according to the invention, 

10 the software, simulates the response of the detector D 
and can therefore produce a spectrum approximately 
identical to this total spectrum. 

The simulation of the response of detector D is 
purely digital and is based on the software that is 

15 stored in the memory 2 6 and in which a number of items 
are input, namely (a) radioactive emission spectra 
representative of some radioelements or their mixes and 
that are also memorized in the memory 26, (b) detection 
characteristics in the form of coefficients and data 

2 0 modelling the thickness through which y radiation 

passes and therefore representing the attenuation, (c) 
operating characteristics of received y radiation, in 
particular representing the detector aperture angle , 
the detected energy bands and amplification of the 
25 electronics and (d) a mathematical motor, for individual 
reproduction of y radiati on emitted for the chosen 
radioelements or the chosen mixes of radioelements. 

The generation of representative radiation counts 
is simulated using a Monte Carlo method using random 

3 0 numbers. 

Thus, simulated responses of the detector D can be 
obtained in order to build up regression straight lines 
similar to those used in prior art but without making a 
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genuine measurement. The detector is simply calibrated 
with an arbitrary rod in the batch to be inspected, for 
which the real composition was analyzed in advance. 
The response of the detector D obtained with this rod 
5 can be used to calculate a fictitious mass that will be 
useful in all subsequent calculations. 

The following description provides additional 
information about the simulation. 

With the software, a count number is calculated 

10 and an attempt is made to reproduce the count number 
for a real pellet. The rod used for the calibration is 
used to "calibrate" the software. The isotopic 

composition and percent of each radioelement in each 
pellet in this rod are known. When the simulation is 

15 done,, the activity of this isotope is calculated for 
each isotope i contributing to the spectrum, taking 
account of all energies j and all attenuations k 
between the scintillator and a fictitious pellet. 

. The mathematical motor makes a gaussian 

20 distribution of the energy as - a function of the 
resolution. 

The result is a simulated spectrum for each 
isotope i considered. A total spectrum for all 

isotopes is obtained by summation. Therefore, the 
25 response of the detector is simulated. The next step 
(or the previous step in another particular embodiment) 
is to use the calibration rod. A true count is 
obtained for each pellet in the calibration rod, using 
detector D. 

30 The simulation is started. The parameters for the 

pellet thus inspected are input into the computer. The 
computer calculates a spectrum for an area of interest 
(in other words one or several energy bands in which we 
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are interested) . If the count obtained by simulation 
is too high, or if it is not high enough, the 
fictitious mass is corrected until the count obtained 
is the same as the count obtained with the actually 
5 inspected pellet (the fictitious mass being the mass of 
the pellet for which it is required to simulate a 
response by the detector, using the amplitude of the 
spectrum rather than its shape) . When the count 
obtained is equivalent, in other words when the correct 

10 fictitious mass is obtained, this value will be saved 
and will be used to calculate all other points on 
regression straight lines. 

Many other detectors could be used instead of the 
annular sodium iodide scintillator detector. For 

15 example, a plane Nal scintillator detector or a GeLi 
scintillator detector could be used. 

The invention can also be used to simulate and 
inspect other fuel elements than the rods mentioned 
above for which the pellets are usually cylindrical. 

20 For example, plate shaped elements containing non- 
cylindrical pellets could be simulated and inspected. 
Furthermore, this invention is not limited to 
simulation and inspection of nuclear fuel elements. It 
can be used to simulate and inspect many other 

25 radioactive objects, for example receptacles made in 
series production and containing a radioactive 
material . 

The following description of- calculation loops 
that could be used to simulate the response of the 
3 0 detector D mentioned above (figures 1 and 2) to a 
nuclear fuel pellet 5 formed by a uranium matrix 
containing several plutonium isotopes (the only 
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isotopes considered in the following calculations), is 
given for guidance only and is in no way restrictive: 

1) Start loop for isotope (i) 



Activity (i) = x M x ti x %isotope (i) x A. (i 



where: Av = Avogadro ' s number = 6.022xl0 23 
A (i) = atomic mass of the isotope (i) 
M = fictitious mass of the pellet 

ti = isotopic content = percentage of Pu in the 
10 matrix 

%isotope (i) = percentage of isotope (i) contained 
in the matrix 

= disintegration constant for isotope (i) 
1.1) Start loop for energy (j) 
15 Resolution^, = Energy (j) x a 

Initial flux (i ,j> = Activity (i) x %emission (i< j} x r|/s 
where: T| = geometric efficiency of the detector 

a = resolution percentage independent of the 
energy, defined experimentally 
20 resolution ( j) = resolution of the detector for 

energy ( j ) 

s = surface area of the detector scintillator 
1.1.1) Start loop for the attenuator element (k) 



1 a(i.j.k) 



6a +6a + Z (k) x6a 

PE(i:j;k> PP (i; j; k ) C (ijk) ^ 



xn (k) x Av/A (k) 



25 where: aa = effective photoelectric absorption cross- 
section PE(ij )k ) ' 

aa = Compton effective absorption cross-section 

C (i,j,k) 
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Ga = Effective cross-section for production of 
pairs' PP (i ,j, k > 

Z (k , = atomic number of the attenuator element 
(k) 

5 p (k , = density of the attenuator element (k) 

A(fc, = atomic mass of the attenuator element (k) 
Attenuation coef f icient (i/ jfk j = exp (-y a(i(j(k) x x (k) ) 
where: X (k) = thickness of the attenuator element (k) 
Final flux (i , j#k) = Initial f lux (if j, k -u x Attenuation 
10 coefficient^) 

1.1.2) End loop for the attenuator element (k) 
with initial f lux (i< j #k _ 1} for k = 1 = initial flux (i/j , 
Nal scintillator: 

M.« = feC + Z (NaI > x<6a->]x - x ft(Na0 

15 Absorption (i<j) {NaI , = final flux (i/j>k) x \jl (Nal) x X (Na i) 

a (i.j> 

^ (Nal) 

where: = effective photoelectric absorption 

PE 0) 

cross-section for Nal. 
Z<NaU = average atomic number of Nal 
oa (Nai) _ compton effective absorption cross- 
20 section for Nal C ( j) 

A(Nai) = average atomic mass of Nal 
P(Nai) = average density of Nal 
X<Nan = Nal scintillator thickness 

1.1.3) Start loop for the procedure to draw values for 
25 the photoelectric and Compton absorption. 

N = 1200 values vn are calculated by drawing random 
numbers according to a Gaussian distribution centered 
on an average value equal to Energy { j, and with a 
standard deviation equal to Resolution ( j } . 
30 Vn = vn x Absorption^ j >NaI 
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Vn- gives absorption spectrum (if j,NaI 

n=l 

1.1.4) End loop to draw values for the photoelectric 
and Compton absorption 
Compton dome : 

5 Energy (j .) = hv (j) / (l+2hv (j) /0 . 511) 
Resolution^ j . } = Energy x a 
Dome scattering (j . } = 

Av 

ffijtf (Nal) x Z (NaIJ x Final flux (i#jfk) x x p (N ai) x 

^(Nal) 

X(Nal) 

10 where: h = Planck's constant 

v ( j) = frequency corresponding to the energy 
denoted Energy ( j), the Compton dome energy 
specific to Nal and denoted Energy being the 
scattered photon energy that is less than 

15 Energy { j) 

odif 

= effective Compton scattering cross- 

c 0'> 

section. 

1.1.5. ) Start loop to draw values for Compton 
scattering . 

20 N=1200 values vn are calculated by drawing 

random numbers using a Gaussian distribution 
centered on an average value equal to 
Energy(j') and with a standard deviation equal 
to Resolution^-) 

25 Vn = vn x Dome scattering (i/ j . ,NaI 

N 

^ vn gives scattering spectrum (i( j } NaI 

1.1.6. ) End loop to draw values for Compton scattering. 
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Compton front (the energy of the Compton front being 
denoted Energy ( j..> and less than Energy (j) ) 

Energy (j ..)=H (2xhv (j) /0 . 511) /[1+ (2xhv (j) /0 . 511) ] r>xhv (j) 
Resolution (j . . ) = Energy (j . . } xa 
Front scattering ( j . . } 

odif 

(NaI)xZ( N ai, x final fluX( i#j(]c , x — xp (NaI) x X (NaI , 

where: adif (Nai) = effective Compton front scattering 
C ( j..) cross-section . 

1.1.7) Start loop to draw values for the Compton front. 

N=1200 values vn are calculated by drawing 
random numbers using a Gaussian distribution 
centered on an average value equal to 
Energy (j-) and with a standard deviation equal 
to Resolution ( j . . ) 

Vn = vn x Front scattering { j . , Nal 



N 



^ vn gives scattering spectrum (i/ j . . ,NaI 

1.1.8) End loop to draw values for Compton front 
scattering. 

Calculation of the Compton background (the energy of 

the Compton background being denoted Energy^..., and 

being less than Energy (j) ). 

Energy (j . . .) = hv< j} / ( l+2hv (j) /0 . 511 ) 

Resolut ion ( j . . . ) : fixed by experience 

Background scattering ( j . . . j 

odif Av 

= (NaI)xZ (Na i, x final flux (i , j/k) x xp {NaI) x X (NaI) 

C 0 ) A (NaI) 

odif 

where: m<i = effective Compton background scattering 
C(j ) 

cross-section. 
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1.1.9) Start loop to draw values for Compton background 
scattering . 

N=12 00 values vn are calculated by drawing 
random numbers using a Gaussian distribution 
centered on an average value equal to 
Energy (j...) and with a standard deviation equal 
to Resolution (j ... ) 

Vn = vn x Background scattering (i< j . . . , , isjal 

N 

^ vn gives scattering spectrum (i/ j . . . } NaI 

n=l 

1.1.10) End loop to draw values for Compton background 
scattering . 

1.2) End loop for energy (j) 

2) End loop for the isotope (i) 

An area of interest for an energy band of interest 
is then chosen (for example from 75 keV to 100 keV) and 

N 

the number of pulses is calculated S = (Vn for 

i.j n=l 

absorption spectrum (i# jjNal) in this energy band is then 
calculated. Therefore we neglected scattering but it 
could be taken into account in another particular 
embodiment . 

T m , which is the percentage of global idle time 
specific to the count system 7 of the detector D is 
then determined experimentally, and Sxt m is calculated. 

All calculations can then be repeated for other 
areas of interest, in other words different energy 
bands . 
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